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SPECinCATION 
HELD 

[0001] The present invention relates to tlic generation of oil and gas exploration and 
production data attributes. 



BACKGROUND 

[0002] The present application claims priority to co-pending U.S. Provisional Paleni 
Application Serial No 60/443,351 filed on January 29, 2003. 

[0003] The goal of hydrocarbon exploration is to find porous and permeable geologic 
10 deposits containing high pore-space saturations of hydrocarbons, under sufficient 

pressure to allow some mode of commercial production. In pursuit of this goal, 
companies, countries and individuals collect and process many types of geophysical 
and geological data. The data is often analyzed to find anomalous zones thai can 
reasonably be attributed to the presence of hydrocarbons. 

15 [0004] The usage of 2D and 3D seismic data anomalies has been a standard practice in the 
petroleum industry since the 1960s. Other geologic and geophysical data anomalies 
have been tried, sometimes successfully, for over a century. These include various 
gravimetric, electromagnetic, chemical, biological and speculative methods. 

[0005] Tlic usage of anomalies for oil and gas detection has been plagued by several 
20 problems. First, most remote sensing anomalies (e.g., a 3D seismic amplitude 

anomaly) cannot be directly tied to a rock property that could be measured in the 
laboratory or using well logs. Much effort is expended attempting to lie observed 
anomalies to known rock responses by modeling the expected attribute response or 
otherwise correlating with a known producing reservoir. This work is often based on 
25 the experience of the practitioner. 

[0006] A second problem is that the anomalies themselves are often evaluated or lied lo 
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response models in a qualitative manner. With qualitative assessment as the basis, 
quantitative, objective and reproducible error analysis has not been possible, 

A third problem is that a basic physical property at work in hydrocarbon reservoirs is 
that both oil and gas are less dense than water. This generally causes oil and gas to 
accumulate up^structure in the pore-space of potential reservoir rocks. The higher 
water saturations are found, generally, down-structure. This separation of saturations 
is driven by gravity. When such a separation of fluid types occurs, flat interfaces, in 
depth, are expected to form. 

This separation causes numerous possible classes of data attribute response. First, 
the hydrocarbon reservoir will have one response for each hydrocarbon type. The 
water-saturated part of the reservoir may have a second data response and the 
interfacial area a third type of attribute data response. In tliis sequence of responses, 
neither the water saturated reservoir response nor the single component hydrocarbon 
saturated reservoir response are expected to vary with structural position. The 
transition from one type of fluid saturation to another, having a different density, is 
expected to occur at a single depth or seismic txavcl time within the attribute data set. 
This change from saturation state to another at a given depth and location should be 
detectable using a quantitative tool. 

[0009] The present invention is designed for the detection, quantification and evaluation of 
20 the depth and location of interface between lighter and heavier saturating fluids as 

exhibited in a data attribute dataset. It is designed to quantify the change in a data 
attribute in the up structure direction from an interpreted water reservoir to the 
hydrocarbon reservoir part of the dataset. The invention has been designed to 
overcome these known problems in the art, 

25 

SUMMARY 

[000 lOJ The method for deriving a GrAZ seismic attribute file entails inputting horizon file 

data and, then, obtaining the gradient of the horizon file thereby producing a horizon 
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vector file. The method next involves inputting attribute file data, indexing from the 
attribute file data at conesponding geographic locations of the horizon file, forming 
an attribute file, and obtaining the gradient of the attribute file thereby producing an 
attribute vector file. The method ends by performing a compilation of the horizon 
5 vector file and the attribute vector to ascertain if changes are in a direction towards a 

surface datum for a narrow time and depth range are detected and measured. 



BRIEF DESCRimON OF TlIE DRAWINGS 

[00011] A specific embodiment of the method will be described by way of example with 
10 reference to the accompanying drawings, in which: 

[000121 FIG 1 is a diagram of the overall method; 

[00013] FIG 2 depicts a geographic depiction of the horizon file within a large geographic 
area having structural contours; 

[00014] FIG 3 depicts the area of the attribute file containing attxibute contours; 

15 100015] FIG 4 (lepicis the hydrocarbon reservoir interface with a water reservoir; and 

[00016] FIG 5 depicts the interface as a plurality of hydrocarbon reservoir types and water 
reservoirs separated by interfaces. 

[00017] The present method is detailed below with reference to the listed Figures. 



20 DER AILED DESCRin^lON OF THE PREFERRED EMBODIMENTS 

[00018] Before explaining the present method in detail, it is to be understood that the method is 
not limited to the particular embodiments and that it can be practiced or carried out 
in various ways. 

[00019] The present method was conceived to detect the changes in a data attribute response 
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in a data set associated with the change from a water reservoir to a hydrocarbon 
reservoir in a given area for a given geologic unit. In addition, the method was 
conceived to operate in high noise, low signal to noise environments, where the data 
quality is poor- The method was designed to operate on subtle hydrocarbon 
5 indicators. The method was also designed to characterize fully the isotropy of the 

water reservoir and of the hydrocarbon reservoir away from the interface between 
these reservoirs. Finally, tlic method was also conceived to determine the errors and 
uncertainties in all measurements and data attribute results relative to a given 
hydrocarbon reservoir and the conesponding water reservoir and the intcrfacial 
10 region. 

[00020] The method addresses the case of multiple hydrocarbon zones, such as, gas qvcr oil 
over water. It is designed to specifically test a model wherein gas is less dense than 
oil and oil is less dense than water, and data responses vary by structural position, 
but transition in a narrow range of depths or two-way seismic times. In a preferred 
IS embodiment, the narrow range would not exceed 5% of the total seismic two-way 

time or depth range contained witliin the horizon file. 

[00021] Within a portion of the data attribute datasel, characterizing a single saturating fluid, 
this method is designed to verify that the data attribute response is invariant wiih 
respect to structural position. The quantification of this invariancc and the associated 
20 uncertainty allow baselines to be established with respect to which the significance 

of an interfacial signal can be assessed. 

This method is also designed to quantify responses and quantify response 
uncertainties in a manner that can be consistently defined, reported and replicated by 
others. Quantification and replication make the output of this method suitable for 
quantitative comparison with rock physics analyses, petro-physical analyses, 
response modeling and geologic analyses (eg., fit to structure analysis). This 
quantification allows this method to be used in those cases where the hydrocarbon 
response is subtle and/or the signal to noise level is low. This combination of 
capabilities represents an advance over current methods. 
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[00023] 



[00024] 



[000251 



[00026] 



15 



20 [00027] 



25 [00028] 



The preferred embodiment h a method for computing a new data attribute for an area 
of interest An area of interest means within the context of this patent application, 
cither a geographic area for a hydrocarbon reservoir, an associated water reservoir, 
contiguous combinations of these or combinations of these with other reservoirs. 

The preferred embodiment also relates to a method for deriving a GrAZ seismic 
attribute file. The method begins by inputting a horizon file and inputting an 
attribute file. 

The attribute file and the hori;con file must contain data elements at corresponding 
geographic coordinates. The geographic coordinates can be an X-Y prospect 
coordinate system, X-Y field development system^ latitude and longitude, internal 
3D seismic survey coordinates, and combinations thereof. The data must be 
sufficiently continuous to allow the computation of a first derivative in each 
coordinate direction in each file, neglecting geographic edge effects. 

The next step involves obtaining the gradient of the horizon file thereby producing a 
horizon vector file, having components, dHl and dH2 at each location G. The first 
component, dHl^ is the partial derivative of the horizon depth or seismic two-way 
time in the first coordinate direction at G. The second component, dH2, is the partial 
derivative of the horizon depth or seismic two-way time in the second coordinate 
direction at G- 

The next step involves obtaining the gradient of the attribute file thereby producing 
an attribute vector file, having components, dAtl and dAt2 at each location G. The 
first componcm, dAil, is the partial derivative of the attribute value in the first 
coordinate direction at G. The second component, dAt2, is the partial derivative of 
the attribute value in the second coordinate direction at G. 

The final step of the method involves performing a compilation of the horizon vector 
file and the attribute vector file to form a combined attribute file, GrAZ. The 
combined attribute file is studied to ascertain if observed changes in the attribute file 
are in a direction towards a surface datum for a narrow range of time and depth. If 
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any components of either the horizon vector file or the attribute vector file do not 
exist or can not be computed at a location G, then no member of the combined 
attribute file, GrAZ exists at G or is assigned to the location G. 

[00029] After the compilation is performed, horizon binning, plotting or other analysis 
5 techniques may be performed to analyze the combined attribute file, GrAZ. 

[00030] The method can further include the step of using dot product matlicmatics to perform 
the compilation- The dot product mathematics is a summation at each geographic 
location G of the product of corresponding elements of the horizon vector file and 
the attribute vector file at each geographic location G. 

10 [00031] For example, if at a location G, the components of the horizon vector file arc dHl 
and dH2 and if at a location G, the components of the attribute vector file are dAtl 
and dAt2, then the dot product is the sum of dHl multiplied by dAtl with dH2 
multiplied by dAt2. If any of the quantities dH1» dB2» dAtl or dAt2 do not exist or 
cannot be computed as finite real numbers at a location G, then a dot product is not 

15 performed and no element of the combined attribute file, GrAZ, exists at G or is 

assigned to the location G. 

[00032] The method contemplates that the horizon file is a lime horizon file made of a set of 
two-way seismic time values depicting the seismic travel time from the datum to the 
horizon of interest and back to a datum. Tn addition, the horizon file can be a depth 
20 horizon file made of a set of values that depict the depth from a datum to the horizon 

of interest within the Earth. 

[00033] The attribute file in the method can one member of the following: 

a. a set of compiled seismic reflection data processed using a defined attribute 
generating algorithm^ and extracted for a horizon of interest; 

25 b, a set of compiled seismic reflection data proces.sed using a defined attribute 

generating algorithm in conjunction with a horizon of interest; 
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c. a set of compiled seismic velocity data processed usiug a defined attribute 
generating algorithm and extracted for a horizon of interest; 

d. a set of compiled seismic velocity data processed using a defined attribute 
generating algorithm in conjunction with a horizon of interest; 

5 e. a set of geophysical gravity data extracted for a horizon of interest; 

f. a set of geophysical gravity data compiled for a horizon of interest; 

g. a set of geophysical gravity data collected for a horizon of interest; 

h. a set of geophysical remote sensing data extracted for a horizon of interest; 

i. a set of geophysical remote sensing data compiled for a horizon of interest; 
10 j\ a set of geophysical gravity data collected for a horizon of interest; 

k. a set of compiled geologic measurements for a horizon of interest; 

I. a set of collected geologic measurements for a horizon of interest; 

mi. a set of petro-physical measurements for a horizon of interest; 

n. a set of compiled or collected engineering data for a horizon of interest; and 

15 o. combinations thereof. 

[00034] The method is better understood with reference to the Figures. 

[00035] FIG 1 is an overall schematic of the method for deriving a seismic attribute file. 

After the initial. step of inputting horizon file data, (Step 110), the gradient of the 
horizon file data is obtained thereby producing a horizon vector file (Step 14-0). The 
20 other initial step of inputting attribute file data (Step 120) is followed by indexing 

from the attribute file data at corresponding geographic locations of the horizon file 
(Step 130). The product of this step is an attribute file. The gradient of the attribute 
file data is then obtained producing an attribute vector file (Step 150). A compilation 
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of the horizon vector file and the attribute vector file is peiformcd (Step 160). 
Finally, from the compilation, it can be ascertained if measurable changes are in a 
direction towards a surface datum for a specific time and depth (Step 170). The 
specific time and depth can then be detected and measured. 

5 [00036] FIG 2 depicts a geographic depiction of the horizon file within a large geographic area 
19 having structural contours 21 and 22. A geographic location G is depicted by 
element 14. At geographic location G, a gradient vector is derived from a horizon file 
for the large geographic area 1.9 and is depicted as element 24. The horizon gradient 
vector 24 has two components, dHl is element 25 and dH2 is clement 26. 

in [00037] The first component, dHl, element 25, is the partial derivative of the horizon depth or 
seismic two-way travel time in the first coordinate direclion at geographic location G. 
The second component, dH2, clement 26, is the partial derivative of the horizon depth 
or seismic two-way travel time in the second coordinate direction at geographic 
location G. 

15 [00038] Element 10 depicts the first coordinate direction and clement 11 depicts the second 
coordinate direction- It is assumed that the coordinate system is orthogonal in a 
preferred embodiment. 

[00039] Tlie axis set illustrates the two coordinate directions and applies to both the horizon file 
and the attribute file. 

2[) [00040] Further, FIG 3 depicis the area 19 of the attribute file containing attribute contours 30, 
31 and 32. The geographic location, G an attribute gradient vector 34 is derived from 
an attribute file for the large geographic area 19- 

[00041] The attribute gradient vector 34 has two components, dAtl shown as element 35 and 
dAt2 shown as element 36- Tlie first component 35, dAil, is the partial derivative of 
25 the attribute file in the Cist coordinate direction at G, shown in FIG 2, The second 

component 36, dAt2, is the partial derivative of the attribute file in the second 
coordinate direction at G, shown in FIG 2. 
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[00042] As an example of compilation, the dot product at geographic loccition G is formed by 
adding (dHl multiplied by dAtl) to (dH2 multiplied by dAt2). If any of the quantities 
dW:U dH2, dAtl or dAt2 do not exist or can not be computed as finite real numbers at a 
geographic location G, tlien no dot product is performed and no element of the 
5 combined aitribule file, gradient attribute depth (which phrase is referred to herein as 

'"GrAZ"), exists or is assigned to the geographic location G- This method of 
computation is a standard mathematical dot product of two vectors at a geographic G. 
The set of all such dot produas successfully computed over the area 19 is 
geographically indexed to form the combined allribule file, GrAZ. 

10 [00043] FIG 4 depicLs the geometry of a hydrocarbon reservoir 38, a water reservoir 40 and the 
interface 42 between the hydrocarbon reservoir 3S and the water reservoir 40. Jt is 
also possible that water reservoir 40 could be a hydrocarbon reservoir and if that is 
the case then hydrocarbon reservoir 38 can be a lighter hydrocarbon reservoir as 
compared with the reservoir 40. 

15 [00044] The datum surface 44 is the datum from which all depths and times arc measured. 

Hlemcnt 50 indicates the direction of increasing depth and/or increasing seismic two- 
way travel lime as measured from the datum surface 44. 

f00045] In FIG 4, a lop surface 46 gives the bounds the hydrocarbon/water reservoirs 38 and 
40. The bottom surface 48 gives the lower bounds of the hydrocarbon/water 
20 reservoirs 38 and 40 respectively. 

[000461 Top surface 46 or bottom surface 48 may serve as a typical horizon of interest. The 
horizon of interest may be any surface defined to be a surface expressible as a one- 
to-one function of either (or both) Ihe lop surface 46 or the bottom surface 48. 

[00047] Tlie geologic region of proposed high water saturation and minimal or no productive 
25 hydrocarbon saturation 40 is shown in FIG 4, The interface 42 divides the two 

geological regions, 38 and 40. The interface 42 is normally expected to be flat with 
respect to the depth or seismic two-way travel time axis. 

[00048] FIG 5 depicts the interface as a plurality of hydrocarbon reservoir types and water 
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reservoirs separated by interfaces. The interface is shown with respect to a depth or 
time axis 50 which indicates the direction of increasing depth and/or increasing 
seismic two way travel time from a surface datum 44. 

[00049] The bottom layer 52 is predominately a water reservoir containing a high water 
saturation and minimal or no productive hydrocarbon saturation. The adjacent layer 
54 is a hydrocarbon reservoir containing a high hydrocarbon saturation of a type that 
IS less dense than the saturating fluid(s) in the bottom layer 52. 

[00050] The middle layer 56 is anotlier hydrocarbon reservoir containing a high hydrocarbon 
saturation of a type that is less dense than the saturating fluid(s) in the adjacent layer 
54. 



[00051] The top layer 58 is a hydrocarbon reservoir containing a high hydrocarbon saturation 
of a type that is less dense than the saturating fluid(s) in the middle layer 56. 



15 



20 



2? 



[00052] It should be noted that FIG 5 is only one embodiment that this method can analyze, 
and various reservoirs can form specific layers of the form shown as 52, 54, 56 and 
58. 

[00053] In discussing this metliod, several terms require consistent definition. The water 
reservoir is a geologic rock formation having both porosity and permeability and 
saturated primarily by water. The water formation may contain a partial 
hydrocarbon saturation, but at a sufficiently low level so as to preclude economic 
development. 

[00054] Similarly, the hydrocarbon reservoir is a geologic rock fonnation having both 
porosity and permeability and saturated in most cases by a combination of water and 
hydrocarbons. The saturation of hydrocarbons must be sufficiently high so as to 
allow economic development. If the saturation of hydrocarbons does not allow the 
production of hydrocarbons and associated water in quantities that arc commercial, 
the reservoir would not be called a hydrocarbon reservoir. Typically, the 
hydrocarbon reservoir is found up-structurc of the water reservoir, which is located 
down-structure. In this discussion, up-siructure refers to shallower depths from the 
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surface within the earth. Down-structure refers to deeper depths within the earth. In 
the case of seismic travel times, deeper depths correspond to larger absolute value 
seismic travel times and shallower depths to smaller absolute value travel times. 

[00055] Both deptlis and seismic travel times are typically measured from a specified datum, 
5 The datum is a specified surface to wliich measurements arc referenced. For 

example, in offshore exploration and production, the datum is typically taken to be 
mean sea level. Depths or seismic times are then referenced to mean sea level as the 
datum. 

r00056] While this method has been described with emphasis on the preferred embodiments, 
JO it should be understood that within the scope of the appended claims^ the method 

miglit be practiced other than as specifically described herein. 
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